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ABSTRACT 

We estimate the contributions to the cosmic microwave background radiation 

O 
O 



(CMBR) power spectrum from the static and kinematic Sunyaev-Zel'dovich (SZ) 
effects, and from the moving cluster of galaxies (MCG) effect. We conclude, 
^ ■ in agreement with other studies, that at sufficiently small scales secondary 

fluctuations caused by clusters provide important contributions to the CMBR. 
At t ^ 3000, these secondary fluctuations become important relative to lensed 
primordial fluctuations. Gravitational lensing at small angular scales has 
been proposed as a way to break the "geometric degeneracy" in determining 
fundamental cosmological parameters. We show that this method requires the 
separation of the static SZ effect, but the kinematic SZ effect and the MCG 
effect are less important. The power spectrum of secondary fluctuations caused 
by clusters of galaxies, if separated from the spectrum of lensed primordial 
fluctuations, might provide an independent constraint on several important 
cosmological parameters. 
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1. Introduction 

The power spectrum of the cosmic microwave radiation (CMBR) carries much 
cosmological information about primordial density fluctuations in the early Universe. 
As photons leave the last scattering surface and travel across the Universe, however, 
these brightness fluctuations are modified by intervening structures, causing secondary 
fluctuations, which we would expect to become more important on small angular scales. 

The power spectrum of the CMBR alone can be used to determine cosmological 
parameters. Recently it has been shown, however, that a geometrical degeneracy effect 
prevents some combinations of cosmological parameters from being disentangled by the 
power spectrum alone ( |Bond, Efstathiou fc Tegmark 1997] ; [Efstathiou fc Bond 19981 



Zaldarriaga, Spergel fc Seljak 1997]) . The primordial density fluctuations and matter 



content determine the positions and magnitudes of the Doppler peaks at the last scattering 
surface. These fluctuations are transferred to apparent angular scales determined by their 
angular diameter distance. As a result, cold dark matter (CDM) models with the same 
primordial density fluctuations, matter content, and angular diameter distance can not 
be distinguished. The models are "effectively degenerate" in the sense that their power 
spectrum is degenerate for parameter determination on intermediate and small scales. 
Although the observed power spectrum also depends on the time variation of the metric, 
via the integrated Sachs- Wolfe effect, this breaks the degeneracy only at large angular 
scales. Unfortunately observations of the power spectrum do not provide strong constraints 
on models at large scales, since this is where the statistics of the data are dominated by the 
cosmic variance due to the fact that we have only one realization of our cosmological model, 
the Universe itself, and we encounter a sampling problem). Therefore we can determine, for 
example, only combinations such as Oo^ 2 and flbh 2 (where Qq and Qb are the z = matter 
and baryon density parameters and h is the dimensionless Hubble parameter). 

It has been noticed that gravitational lensing can break the geometric degeneracy 
at small angular scales, I ^ 2000, in such a way that the cosmological parameters can 
be determined separately ( Metcalf fc Silk 1998| ; |Stompor fc Efstathiou 19981 ). The effect 



of gravitational lensing on the CMBR was studied by several authors (see for example 
Blanchard fc Schneider 1987] ; Metcalf fc Silk 1997) ; |Seljak 1997] ) . Static gravitational lenses 



do not change a smooth CMBR, but the fluctuations get distorted by lensing. As a result, 
power from the acoustic peaks is transferred to small angular scales, conserving the variance 
of the spectrum. The amount of power transferred depends on the cosmological model, 
thus, in principle, we can determine separately Qq, and h. This method makes use of 
the small angular scale part of the power spectrum, where the amplitude of primordial 
fluctuations is declining and secondary fluctuations are becoming more important. The 
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question naturally arises: How do contributions to the power spectrum from secondary 
fluctuations influence parameter determination based on the small scale CMBR power 
spectrum? 

The most important secondary fluctuations are the thermal static and kinematic 
Sunyaev-Zel'dovich (SSZ and KSZ) effects ( (Sunyaev fc Zel'dovich 1980|) , the Rees-Sciama 
(RS) effect QRees fc Sciama 1968| ) , the moving cluster of galaxies (MCG) effect ( pirkinshaw 



|fc Gull 1983| ; |Gurvits fc Mitrofanov 1986| ; |Pyne fc Birkinshaw 1993]) , point sources 



( [Toffolatti et al. 1999|) , and, if the Universe was re-ionized at some early stage, the 



Ostriker-Vishniac effect ( Pstriker fc Vishniac 1986| ; |Vishniac 1987| ). In this paper we 



concentrate on secondary effects caused by clusters of galaxies. Since detailed reviews are 
available on the SZ and the MCG effects ( [Rephaeli 1995| ; [Birkinshaw 1998| ), here we just 



summarize their major features. The thermal SZ effect is a change in the CMBR via inverse 
Compton scattering by electrons in the hot atmosphere of an intervening cluster of galaxies. 
We use the terms kinematic or static SZ effect depending on whether or not the intracluster 
gas possesses bulk motion. To date only the static thermal SZ effect has been detected 
( |Birkinshaw"T998|) . 

The MCG effect is a special type of RS effect, due to the time- varying gravitational 
field of a cluster of galaxies as it moves relative to the rest frame of the CMBR. Unlike the 
original RS effect, the MCG effect is not caused by intrinsic variation of the gravitational 
field, so that in the rest frame of the cluster, the photons fall into and climb out of the same 
gravitational field. However, in the rest frame of the cluster the CMBR is not isotropic, but 
has a dipole pattern, being brighter in the direction of the cluster peculiar velocity vector. 
Photons passing the cluster are deflected towards its center. Thus in the direction of the 
cluster peculiar velocity vector (ahead of the cluster) one can see a cooler part of the dipole 
pattern. Towards the tail of the cluster, one can see a brighter part of the dipole (ahead of 
the cluster) . When transferring back to the rest frame of the CMBR, we transfer the dipole 
out, but the fluctuations remain, showing a bipolar pattern of positive and negative peaks. 
At cluster center there is no deflection, thus there is no effect. The amplitude of the MCG 
effect is proportional to the product of the gravitational deflection angle and the peculiar 
velocity of the cluster. 

The most important characteristics of the SSZ, KSZ, and MCG effects in the context 
of cosmology is that their amplitudes do not depend on the redshift of the clusters causing 
the effect. Using thermodynamic temperature units, their maximum amplitude are about 
500 /iK, 20 fiK, and 10 /iK, respectively. The SSZ and KSZ effects have the same spatial 
dependence as the line of sight optical depth, the MCG effect has a unique bipolar pattern. 
Assuming a King approximation for the total mass and an isothermal beta model for the 
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intracluster gas, the full width at half maximum (FWHM) of the SSZ and KSZ effects 
~ (2 — 4) r c , where r c is the core radius, depending on (3, which is typically between 2/3 
and 1. The MCG effect has a much larger spatial extent, with FWHM for each part of the 
bipolar distribution « 10 r c . 

The spectra of the effects are also important: the SSZ effect has a unique spectrum 
which changes sign from negative to positive at about 218 GHz. The KSZ and MCG 
effects have the same frequency dependence as the primordial fluctuations. The most 
important difference between the SZ effect and the MCG effect is that the SZ effect is 
caused by intracluster gas, the MCG effect is caused by gravitational lensing by the total 
mass regardless the physical nature of that mass. Therefore the SZ effect only arises from 
clusters with intracluster gas. Clusters can produce significant MCG effects even if devoid 
of intracluster gas. 

The effects of clusters of galaxies on the CMBR in a given cosmology have been 
a subject of intensive research since the late 1980s. There are several different ways of 
extracting information from these effects. Source counts of the SSZ effect were estimated 
by using the Press-Schechter mass function (PSMF) and scaling relations ( |Cole fc Kaiser 



158$ IMarkevitch et al. l992| 1994; |Makino fc Suto M: [Bartlett fc Silk 1994 De Luca 



Desert & Puget 1995, Colafrancesco et al. 1994j ; 1997, [3uto et al. 1999|) . The importance 



of the SSZ effect was demonstrated and it was shown that thousands of detections are 
expected with the next generation of satellites. Contributions to the CMBR from the RS 
and the KSZ effects were derived by Tuluie, Laguna & Anninos (1996) and Seljak (1996) 
for CDM models with zero cosmological constant. Tuluie et al. used N-body simulations 
and a ray-tracing technique, Seljak used N-body simulations and second order perturbation 
theory. Contributions from the SSZ and KSZ effects originating from large scale mass 
concentrations (superclusters) were studied by Persi et al. (1995). Bersanelli et al. (1996), 
in their extensive study of the CMBR for the Planck mission, estimated the contribution 
to the power spectrum from the SSZ and KSZ effects. Aghanim et al. (1998) estimated 
the effects of the KSZ and MCG effects on the CMBR including their contributions to 
the CMBR power spectrum. Aghanim et al. simulated 12.5° x 12.5° maps with pixel 
size of 1.5' x 1.5'. They used the PSMF normalized to X-ray data (assuming an X-ray 
luminosity- mass relation). The total mass was assumed to have a Navarro- Frenk- White 
profile ( |Navarro, Frenk fc White 1997| ), and the intracluster gas was assumed to follow an 
isothermal beta model distribution. The time evolution of the electron temperature and 
the core radius were assumed to follow models of Bartlett & Silk (1994), which are based 
on self-similar models of Kaiser (1986). According to Aghanim et al. (1998)'s results, the 
KSZ effect is many orders of magnitude stronger than the primordial CMBR on small 
angular scales, and therefore the effect would prevent the use of the power spectrum to 
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break the geometric degeneracy. Atrio-Barandela & Mucket (1998) estimated the power 
spectra of the SSZ effect in a standard dark matter dominated model with different lower 
mass cut-offs. Contributions to the power from the Ostriker-Vishniac effect in CDM models 
were estimated by Jaffe & Kamionkowski (1998). 

In this paper we estimate the contributions to the CMBR power spectrum from the 
SSZ, KSZ, and MCG effects on small angular scales adopting cold dark matter dominated 
models. In our models we assumed scale invariant primordial fluctuations with a processed 
spectrum having a power law form on cluster scales with a power law index of rip = —1.4 
(Bahcall & Fan 1998). This maybe used as a first approximation as long as contributions 
from very low and/or very high mass clusters are small (cf. our discussion about mass cut 
offs at section []). We use three representative models in our study: Model 1, open CDM 
(OCDM) model: a low density open model with Q = 0.2, A = 0, cr 8 = 1.2; Model 2, flat 
lambda CDM (ACDM) model: a low density flat model with Q = 0.2, A = 0.8, cx 8 = 1.35; 
Model 3, standard CDM (SCDM) model: a flat model with f2 = 1, A = 0, and cr 8 = 0.65. 

In Section ^| we outline our method of estimating the power spectra of secondary 
fluctuations caused by clusters of galaxies and discuss our normalization method for the 
PSMF. Sections |3| and ^ describe how we used the PSMF and the scaling relations to obtain 
masses and other physical parameters of clusters. In section ||] we present the spherical 
harmonic expansion of the SSZ, KSZ and MCG effects, and our method of estimating their 
power spectra. Section ^| describes our simulations to evaluate the integrals over clusters. 
Sections [5] and |8] present our results and discuss the differences from previous work. 



2. Outline of the method 



We used the Press-Schechter Mass Function (PSMF, |Press fc Schechter 1974|) as a 



distribution function for cluster masses. We used n = —1.4 as indicated by observations 
(Bahcall and Fan 1998, and references therein). We used observationally determined 
cluster abundances as a constraint on the PSMF. Where necessary, we altered the model 
parameters resulting from the usual top hat spherical collapse model since that model is 
only an approximation. In our SCDM model we changed only the overall normalization of 
the PSMF by multiplying it by 0.23 (a similar normalization was used by |De Luca et al.| 
1995). This procedure is inconsistent with the interpretation of the PSMF as a probability 



distribution (it does not integrate to unity), but we use results for the SCDM model only 
as a comparison to the other two models. In our Lambda-CDM model we multiplied the 
critical density threshold, 5®, obtained from the spherical collapse model (equation @), by 
1.23 (which is equivalent to changing the o~g normalization) and made no other changes. 
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Our OCDM model needed no adjustments. With these changes all three models agree well 
with the present day (z = 0) observed mass spectrum (Figure [I]; Bahcall & Cen 1993). 
For high masses, the first two models (OCDM, and ACDM) also agree with the observed 
z dependence of the high mass cumulative mass function (Figure ^ Bahcall, Fan & Cen 
1997). The ACDM and the SCDM models agree with CMBR constraints, while the OCDM 
model is rejected by these constraints flLineweaver fc Barbosa 1998|) . As a cautionary note, 
it is useful to keep it in mind that taking all data into account none of these models are 
acceptable. 

We assumed that the total mass distribution follows a truncated King profile ( |Kingj 
For the intra-cluster gas we assumed an isothermal beta model (Cavaliere & 



Fusco-Femiano 1976). Isothermal beta model fits to X-ray images of clusters give f3 ~ 2/3 
( pones fc Forman 1984] ) . Determinations of the (3 parameter based on spectroscopy suggest 
/3 ~ 1 flGirardi et al. 1996 ; Lubin fc Bahcall 1993Q . Numerical simulations imply a range 



for P from 1 to about 1.3 {(3 = 1.05: |Evrard, Metzler fc Navarro 1996|) ; (3 ~ 1.3: Bryan & 



Norman 1998| ; [Frenk et al. 1999] suggest (3 = 1.17). We follow the precepts of [Eke et al. 
1998 and adopt (3 = 1. This choice provides a mass temperature function which is in a 
good agreement with the observed function ( [Horner et al. 1999| ). The fitted X-ray spatial 
distribution is highly dependent on the X-ray structure of the core, and may be expected 
to be less reliable in the outer regions of clusters. The SZ effect is more sensitive to the 
outer regions (the SZ effect is proportional to the electron density as opposed to thermal 
bremsstrahlung, which is proportional to density squared). Choosing the spectroscopically 
derived (3 = 1 gives a smaller SZ effect, and so our choice of (3 should provide a conservative 
estimate of the contribution of the SZ effect to the power spectrum. A slightly larger /3 
would not change our results significantly. Although the beta model describes the inner 
intra-cluster gas well, for more accurate SZ work an improved cluster model, which fits the 
outer regions better, will be needed. 

The other physical parameters were determined using the virial theorem, a spherical 
collapse model, and models of the intra-cluster gas by Colafrancesco & Vittorio (1994). 
We assumed a Maxwellian distribution for the peculiar velocities, v pec , and used results 
of N-body simulations by Gramann et al. (1995) to normalize the distribution. We 
took velocity bias into account, and assumed that the peculiar velocities are isotropically 
distributed. 

We used analytical approximations to calculate the contributions of the SSZ, KSZ and 
MCG effects to the CMBR power spectrum. These contributions are important only at 
small angular scales, where we can neglect the overlap between cluster images and ignore 
the weak cluster-cluster correlation, and therefore we can approximate the resulting power 
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spectrum by summing the contributions from individual clusters (similar methods were 
used by [Cole fc Kaiser 1988| ; [Bartlett fc Silk 1994| ; and |Atrio-Barandela fc Mucket 1998Q . 
We expanded the SSZ, KSZ, and MCG effects as Laplace series (i.e., series of spherical 
harmonics), then determined the individual cluster contributions and summed over the 
clusters (for a detailed description see [Molnar 1998|) . 



Our approximation breaks down at large angular scales, but at these scales primordial 
fluctuations dominate the CMBR, and the cluster contribution is only a minor perturbation, 
so that only a rough indication of the cluster effect is needed. 



3. Total Mass Distribution 

According to the Press-Schechter method, the co-moving number density of clusters of 
total mass M at redshift z (the PSMF) is given by 



dn c (M,z) 2 VLqp c / d In a 



where f2 is the matter density today in units of the critical density, p c = 1.88 x 10 29 gem 3 
is the current critical density of the universe (we adopt a dimensionless Hubble parameter 
hwo = 0.5 in our work); v{M, z) = S^(Qq, z)/a(M), where the present mass variance for a 
power law power spectrum, P{k) oc k n , is 

a — (n + 3)/6, M 8 = 6 x 10 14 f2 h~ 1 M Q is the mass within an R 8 = 8h _1 Mpc sphere, and 
o"8 is the normalization (Lacey & Cole 1993, Press & Schechter 1974). The over-density 
threshold linearly extrapolated to the present may be expressed as ( [Lacey fc Cole 19~93| ; 
Navarro, Frenk fc White 1997|) 



; (12tt) 2 / 3 (1 + 



6° c (n ,z) = { §D(fio,0) 



2tt 



,2/3 



+ 1 



sinh rj—ri > 

0.15(127r) 2 / 3 ^ 0055 ^A(fio, 0)/D A (fio, z) 



where the conformal time for open models is 



Q = 1, A = 

n < i, a = o 

Q < 1, A = 1 - Q 



(3) 



cosh 



^m(^) 



- 1 



(4) 
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For open models with cosmological constant A = the linear growth factor is (Peebles 
1980) 



D(Q , z) = l + ^ + 3(1+ 3 ^ )1/2 ln((l + wf" - 
where w(Qq, z) = Q^-(z) — 1, and the density parameter Q m (z) is 



(5) 



fto(l + z) 



SZo(l + *) + (1 - fio) 

For models with a non-zero cosmological constant the integral can not be done analytically, 
thus we use an approximation (Lahav et al 1991; Caroll, Press & Turner 1992) to obtain 
the equivalent expression to (M, 



D A (Q , z) = (l + z)- l5 -^Uf(zf 7 - n A (z) + 



1 + 



n f (z) 



i + 



n A ( z ) 



70 



(7) 



The density parameters, £lf(z) and Q\(z), for spatially flat Universes with A = 1 — Q , are 



«/(*) = nJl^U (8) 

The normalization of these growth functions is chosen so that at high redshifts they 
approximately match the time variation of density contrast in an Einstein-de Sitter (fl = 1) 
Universe, which is a good approximation to the early Universe whatever its density 
parameter today. The total number of clusters at redshift z (in a redshift interval of dz) is 

, . [Mu P dnJM, z) dV 

N(z)dz = / cV , ' ; — dMdz, 10 

v ; Jm 1ow dM dz ( J 

where Mi ow and M up are the lower and upper mass cut offs for clusters. We used 

M low = 10 13 M for the SSZ and KSZ effects, M low = 10 12 M for the MCG effect, and 

M up = 1 x 10 16 M for all effects. The lower cut off, M [ow , for SZ effects signifies the lowest 

cluster mass for which we expect a well-developed intracluster atmosphere. In the case of 

the MCG effect, Mi ow is the mass limit from which we consider a mass concentration as a 

cluster ("formation"). We found that low mass clusters do not contribute substantially to 

the power spectrum, thus the lower cut off, M tow , has little effect on our results. The upper 

cut off, M up , has no effect on our results (as long as it is large enough): the probability of 
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getting such a large cluster is negligible, so that the contribution from more massive clusters 
is negligible. 



4. Other physical parameters of clusters of galaxies 

We assumed a truncated King profile for the total mass distribution 



3 

2 \ 2 



p(r) = ! Po^ + k) r<pr c 
I r > pr c 

where r c is the core radius, pr c is the cut off, and an isothermal (3 model for the intra-cluster 

gas 



n r 




.M 

r <pr c ; ( 12 ) 
r > pr c 

where n(r) and no are the electron number density at radius r and at the center of the 
cluster (Cavaliere & Fusco-Femiano 1976). Analytical studies and numerical simulations 
show that the gas density profile scales with the total density, and that the gas central 
electron density may be expressed as 

m p (l + X) 

where X = 0.69 is the average Hydrogen mass fraction, f g is the intra-cluster gas mass 
fraction, po, the central mass density, is determined from the total mass by integrating 
equation (|TT|). Little is known about the total mass and redshift dependence of the 
intra-cluster gas from observations. Here, we adopt Colafrancesco & Vittorio (1994) 's model 
which assumes that changes in the intra-cluster gas are driven by entropy variation and/or 
shock compression and heating. According to their model, the gas mass fraction may be 
expressed as 



'•-'"fcp^)'^-' (14) 

where the normalization, f g0 = 0.1, is based on local rich clusters, and we used r\ = 0.5 and 
s — 1, which are consistent with available data. Using the virial radius to express the core 
radius, r c = R v /p, and assuming spherical collapse, we obtain 



-10- 



r c (fi , M, z) = — 

V 

where A c (z) = p v (z)/pb(z) is the overdensity of the cluster relative to the background 
flColafrancesco & Blasi Tggg ). For A = models (our OCDM and SCDM models) the over 
density may be expressed as 



M 



178 



I/O 



lO 15 /i- 1 M y fi A c (fi ,^) 



l + z'' 



(15) 



A c (n, z) = 4vr 2 Q 2 [(Q 2 + 2Q) 1 / 2 - ln(l + Q + (Q 2 + 2Q) 1 ' 2 

where Q = 2(1 - n )/(n (l + z)) (Oukbir & Blanchard 1997). For spatially flat models 
with finite cosmological constant (our ACDM model) we have 



(16) 



A c (tt,z) = 18tt 2 



1 + 0.4093(fy(; 



\ 0.9052 



(17) 



where we used the approximation of Kitayama & Suto (1996). 



Numerical models of cluster formation show that cluster temperature scales with total 
mass. Using the virial theorem and assuming spherical collapse with a recent-formation 
approximation in a standard CDM model, the electron temperature, T e , becomes 



where A c is the density contrast of a spherical top-hat perturbation relative to the 
background density just after virialization (cf. for example Eke, Cole & Frenk 1996). The 
recent-formation approximation, however is valid only for £2 = 1. For our low matter 
density open model (Model 1, OCDM), we use a model which takes into account accretion 
during the evolution of clusters, and leads to the following scaling: 



1+^keV, 



kT e = 2.76 p 



2/3 



2vr 



sinh 7] — rj 



2 / 3 Tip + 3 
+ " 



keV. 



(19) 



This was derived for open models with zero cosmological constant ( Voit fc Donahue 1998| ). 
but since structure formation evolves similarly in a low density model with the same 
matter density and a zero cosmological constant, we use it for our Model 2 (ACDM) as an 
approximation. 



We assumed a Maxwellian for the cluster peculiar velocity distribution, 
expected from a Gaussian initial density field: 



->pea 



as 
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P{Vpeci dVp ec oc v 2 pec exp{-vl ec /2a p (z) 2 }dv pec , (20) 

where u p {z) is the Maxwellian width of the peculiar velocity distribution. The rms peculiar 
velocity from linear theory, smoothed with a top-hat window function of radius R, Wr, is 
given by 



(v 2 )r(z) = H 2 (z) a 2 (z) f 2 (Q , A) cr_i(i?), 
where a(z) is the scale factor, and the moments, <Jj(R), are defined as 



(21) 



a 3 {R) 



2tt 2 



k 2j+2 P(k)W(kR)dk, 



(22) 



where P(k) is the Fourier transform of the power spectrum and equation (plD uses 
the moment j — — 1 (Peebles 1980Q . The velocity factor, f(z) = din 6/ din a, can be 



approximated as (Peebles 1980| ; 1984) 



n f (z) 



A = 
A = 1 - Q 



(23) 



_2(l+z)D A (n ,«) 

The cluster peculiar velocity rms differs from this since we assume that clusters form at the 
peaks of the density distribution, and with this bias may be expressed as 



<«;>*(*) = (v 2 ) R {z) 



(24) 



aKR)v 2 -i(R)- 

(Bardeen et al 1986). Colberg et al. (1998) calculated the velocity bias in a series of CDM 
models using a top-hat filter and processed CDM power spectra. The correction factor has 
a weak dependence on Q : it is about 0.8 for low density and flat CDM models. 



We obtain the Maxwellian width in equation (|20|) from the rms peculiar velocity from 
averaging a Maxwellian: 



f™v 4 exp{-v p ! ec /2a p (z) 2 }dv 
S™v*eM-vlj2v P {z) 2 }dv 



3< 



(25) 



We expressed a p as o p = norm x [H(z)a(z)f(z))]/[H(0)a(0)f(0)]. The normalization at 
z = was determined by using results on the peculiar velocity distribution from numerical 
simulations (|Cramann et al. 1995|) . Thus we obtain the following expression for the 
Maxwellian width of the peculiar velocities, a p (Qo, A, z), with velocity bias for models with 
no cosmological constant (OCDM, SCDM, A = 0): 
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a p (n o ,0,z) = (410kms" 1 )^ 6 (2)(fio(l + ^) + l-^o) 1/2 - 
For our ACDM model (A = 1 - Q ) we obtain 



(26) 



(j p (Qo, 1 — Qq, z) = 410 km s 



-l 



l + z 



x 



(Q (l + zf + l-Q )^ D A (Q ,z) 

5 - 3(1 + z)D A (n , Z y 
5-3L> A (fi ,0) 



x 



(27) 



This normalization is significantly larger than some recent measurements suggest (Bahcall 
& Oh 1996), but it is a good match to others (Gramann 1998). This uncertainty should be 
remembered when interpreting our final results. 



5. Power Spectra of SSZ, KSZ and MCG Effects 

Ignoring the correlation between clusters, the power spectrum becomes 



c 



J dz J 



dM 



dn c (M, z) x dV 
dM 1 dz' 



(28) 



where Gf is the contribution from clusters with total mass M at redshift z, and X denotes 
the SSZ, KSZ or MCG effects. dV/dz is the differential volume element (assumed isotropy) 



dV 

dz 



r(z) 



Attc 



fi (l + zf + (1 - Q - A)(l + zf + A 



'1/2 



(29) 



where the effective distance r(z) is 



r(z) 



A = 



^()^+ffi()-2)(vi+n z-i) " 

-§- o / Z dx[n (i + x) 3 + i - n ]- 1/2 A = i-n 



2c 
^0 



(30) 



(Peebles 1993). In general, the coefficients may be determined by calculating the 
spherical harmonic expansion of the cluster image by averaging out the azimuthal parameter, 
m, 



2£ + 



- 1 m=-i 



(31) 
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Our task is to determine the a^ m coefficients. 

The SSZ and KSZ effects are cylindrically symmetric for spherical clusters, therefore 
we may describe them using only one coordinate, the angular distance from the cluster 
center. We separate the effects into amplitudes and geometrical form factors which carry 
their spatial dependence. The SSZ and KSZ effects in thermodynamic temperature units 
may be expressed as 




x,9) = A s = A° s {x)C(e) 
\ . K x,9)=A K = A° K (x)((e), 
where the central effects for the SSZ and KSZ effects are 



(32) 
(33) 



and 



A° s (x, 9) = [Y (x) + eY x (x) + ® 2 Y 2 (x) + <d 3 Y 3 (x) + Q A Y A (x)} 6 r , 



(34) 



A° K (x,Q) 



rl -r + e(^n + Iyi)} P 2 -[i + ed(x) + e 2 c 2 (x)}pp 1 (a) 



+ [D (x) + QD 1 (x)]/3 2 P 2 (a)\r . 



(35) 



In these expressions Pk is the Legendre polynomial of order of k, x = hv jksTcB is the 
dimensionless frequency, = ^T e / (m e c 2 ) is the dimensionless temperature, a is the angle 
between the cluster's peculiar velocity vector and its position vector, h, u, k B and T C b 
are the Planck constant, frequency, Boltzmann constant, and temperature of the CMBR, 
Tqb = 2.728 ± 0.002 K ( |Fixsen et al. 1996|) , and the lengthy expressions for the spectral 
functions Yi(x), Ci(x) and D^x) may be found in Nozawa, Itoh & Kohyama (1998). 
These functions arise from an expansion of the Boltzmann equation and although they are 
inaccurate for high temperature clusters, their precision is sufficient for our purposes here 
(for a discussion see |Molnar fc Birkinshaw 19UU| and references therein). The optical depth 



through the cluster center for gas model ([12]) is 



r = 2a T n r c I p (3(3/2,p), 



(36) 



and the geometrical form factor is 
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C(9)= (n e /n )dz= (l + ^)3-7 j(J3,p,d c ), (37) 



where the function j(fi,p, d c ) is defined as 



1,(3/3/2, x /p 2 - d*) 

= wa,) ■ (38) 

c? c = 9/9 c holds in the small angle approximation, and the integral, I p (a, q), is 

2 1 (a) 1 (za) 

where T is the gamma function, F(x; y; u>; z) is Gauss' hyper-geometric function, and a 
must be greater than 1/2 flGradshteyn fc Ryzhik 1980|) . The geometrical form factor is 
normalized to one at the cluster center (C(0) = 1). We may expand the SSZ and KSZ effects 
in Legendre series as 



00 IP + 1 

&sz = A° sz £ -1-CiPi, (40) 

where Q are Legendre coefficients of (, and SZ refers to the SSZ or the KSZ effect. We 
determine the Legendre coefficients using a small angle approximation, as 



= 2tt [ Pc ((6)J [(£+-)6}6d6, (41) 
jo I 

where we used the approximation 

P e ~M(t+l)0], (42) 

where J is a Bessel function of the first kind and zero order ( Peebles 1980Q . We can 
convert Legendre coefficients to Laplace coefficients by expressing the Laplace series of such 
a function as 



r?f + 1 IP 4- 1 

£a*,I?(0,yO =EJ^M = E^— ^0*), (43) 

where /x = cos#, and Y™ and are the spherical harmonics and Legendre polynomials. 
Therefore the conversion can be done as 
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Thus the Laplace series of the SSZ and KSZ effects become 

*sz{0) = ALE(^ 1 ) 1/ V/W- (45) 

Using equation ([45]), the contribution of one cluster to the power spectrum of the SSZ 
and KSZ effects becomes 

Gf = ^szfCl (46) 
where A° sz and Q are given by equations fl34]) , (|35"1), and (p]). 
Similarly, the MCG effect may be expressed as 

/ AT" 



V T 



| (x,e,<p) = A M = A 1 Xrffl,<P), (47) 

M 

where £(#, (p) is the geometrical form factor. 9 is the angle of the line of sight relative to the 
center of the cluster. The azimuthal angle, if, is measured in the plane of the sky from the 
direction of the tangential component of the peculiar velocity. The maximum of the MCG 
effect is 



A™ ax = -(v p /c)sma5 max , (48) 

where 5 max is the maximum deflection angle, a is the angle between the cluster's peculiar 
velocity vector, v p , and its position vector, and c is the speed of light in vacuum. For a 
spherically symmetric thin lens the deflection angle is given by 

where d is the impact parameter at the source, and M. is the mass enclosed by a cylindrical 
volume with axis parallel to the line of sight and radius equal to the impact parameter d (cf. 
for example Schneider, Ehlers & Falco 1992). Using the King approximation for the density 
distribution (equation pTH), the total mass in cylindrical coordinates, (R,if),z), becomes 

/•2tt rd c rp RdR 

M=p r 3 c dip dz r , (50) 

Jo J-d c Jo h + R2 + z 2\2 
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where d c = djr c rs 9/9 c in the small angle approximation. A straightforward integration 
and equation (f49|) lead to 



5{d c ,p) 



AGM g{d c ,p) 



c 2 r c g(p,p) d c 



where the function g(x,p) is 



(51) 



g(x,p) = ln(l + x 2 ) + In 



p + y/l+p 2 



P + y/l+p 2 + X 2 

Thus the geometrical form factor in our case becomes 



(52) 



c(a \ 9(dc,p) 



b-mdc 



where b m is the maximum value of the function g(x,p)/x. 



(53) 



The MCG effect depends only on cosy?, therefore we need to determine only the 
m = ±1 terms in the spherical harmonic expansion. Expressing the spherical harmonics by 
associated Legendre polynomials, equation (|5^) expands as 



where we used the identity 



21+1(1-1)! 
4vr (i + 1) 2 ? ^ 



(54) 



(55) 



(£ + m) 

In order to obtain a real function, the imaginary terms must vanish, therefore we must have 



-1 



(56) 



Using orthogonality, expressing the spherical harmonics in terms of associated Legendre 
polynomials and using equations fl56|) and fl5"4"D , we obtain 



in 9{dc,P) n \, a \ 2 

d\l%> (cost)) cos if, 

n*/ d c 



(57) 



where 
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The integral can be performed since g(x,p) and P/ do not depend on ^> giving 



(£-1)! 

4tt ^(£+1)!' 1 



g = -^ r^^p/M, (59) 



from which we find 



C] = + 1/2) # c / <?(M Jx[(£ + -) (9] d0. (60) 

Here we used the small angle approximation for the associated Legendre polynomials: 

Pfo)M (1 + 1/2) J&1 + ±)0\, (61) 

where J\ is the Bessel function of the first kind and order 1 (for a derivation see Appendix). 
Thus the Laplace series of the MCG effect is 

where £j is given by equation ( jBTf ). 

For the power spectrum of the MCG effect, using equation flB^), we obtain 

Cf = ^-(A^flOT. (63) 

The observed effects are calculated by convolving the theoretical fluctuation pattern 
with the telescope's point spread function (PSF). One advantage of using the spherical 
harmonic coefficients is that this convolution is just a multiplication in spherical harmonic 
space. Assuming an axially symmetric PSF, R, its Legendre polynomial expansion may be 
expressed as 

0£+ 1 

R(x ■ x') = ~^RiPi{x ■ x'), (64) 

where the unit vectors, x and x', point to an arbitrary direction (where we want to evaluate 
the expansion) and to the center of the PSF. Assuming a non-axially symmetric effect, its 
spherical harmonic expansion can be written as 
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/(i)=DrW, (65) 

i,m 

where i runs from zero to infinity and m runs from -£ to Using the addition theorem 
for spherical harmonics and their orthogonality, the convolution of these two functions, 
M = R-k f, becomes 



M(x)=J2R e frmx). (66) 

£,m 



6. Simulation of clusters of galaxies 

We used Monte Carlo simulations to generate an ensemble of clusters of galaxies with 
masses sampled from the PSMF (equation [I]) with parameters those of our OCDM, ACDM 
and SCDM models. We obtained the central electron number density and temperature, and 
the cluster core radius from scaling relations (equations [T3|, [T8|, [T9], and |T5] ). 



We choose to sample the PSMF using a rejection method. The magnitude of the 
peculiar velocity may be sampled using an inversion method on the Maxwellian (|20|), and 
yields 

v p = V2a pl -\3/2, RN)-\ (67) 

where 7 _1 (order, x) is the inverse of the incomplete gamma function and a v can be 
determined by using equations (^), (^§) and fl27j). RN is a uniformly distributed random 
number in (0,1). We assumed an isotropic distribution in space for the directions of the 
peculiar velocity vectors, and ignored correlations between cluster peculiar velocities. The 
tangential and radial peculiar velocities are distributed as projections of equation ( |67|) 
accordingly. As an illustration, in Figure |] we show results from one simulation using our 
SCDM model projected on a grid. 

The observational mass function (Figure []]) is specified by M L5 , the mass contained 
within co-moving radius of 1.5 Mpc. To convert Mi. 5 to the virial mass, M v , which we use 
in the PSMF, we assume that the mass profile near 1.5 Mpc can be approximated with 
M V (R) oc R q . We obtain 

w = [ A e (0) \^ 

v ^A c fi m (^)A c (0)(47r/3) Pc (0)(1.5^ 1 Mpc)3 y l L5 ' 1 > 
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or, substituting the numerical values, 



M, 



178 



i 

3-q 



M L 



(69) 



\A c n m (z) 6.98 x W^Mqj 
We used q = 0.64 flCarlberg, Yee fc Ellingson 1997| ) to obtain curves shown in Figure 0. 

The power spectrum for an ensemble of clusters may be determined by summing the 
individual contributions of the simulated clusters (equation |^S|). We binned clusters by 
their apparent core radii, 9 C , then we summed the amplitudes in each bin. The numerical 
evaluation of integral (p8|), in this case, can be performed as 



r*sz 



^" i cl 



2 

SZ)cli 



(70) 



and 



a 



M 



2£ + 



T £(£])< £(Aro 



max\2 
cli 



(71) 



where the index cl runs over clusters whose core radii fall within the i bin 



7. Results 

Our results for the power spectra (more exactly the dimensionless C(£) = £(£ + l)Cg) 
from our Model 1 (OCDM) are shown on Figure [|. Figures |5| and [] show our results for 
Model 2, (ACDM) and Model 3, (SCDM) respectively. As a comparison, in each figure, 
we plot the corresponding primordial CMBR power spectrum (solid line) with CO BE 
normalization including the effects of gravitational lensing calculated by using a new version 
of CMBFAST QZaldarriaga fc Seljak 199"^ ; |Zaldarriaga, Spergel fc Sel.jak 19971 ). On large 



angular scales (up to about £ ~ 10) the cosmic variance dominates (not shown). On small 
angular scales the shape of the power spectra depends on the apparent angular sizes of the 
clusters and the amplitudes of the effects. The apparent angular size depends on how the 
core radius and the angular diameter distance change with redshift, while the amplitude 
is sensitive to the gas content, gas temperature and total mass as a function of redshift. 
Figures f§-|)] demonstrate that for small angular scales {£ ^ 3000) the contribution to the 
power spectrum from the SSZ effect exceeds that of the primordial CMBR in all our models. 
The contributions of the KSZ and MCG effects become important only on small scales, but, 
at those scales, they may dominate over the lensed primordial fluctuations. 
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Due to the early structure formation, there are more clusters at high redshift in our 
OCDM and ACDM than in our SCDM simulations. Therefore the contribution to the 
power spectrum from clusters in a low matter density model is substantially larger than in a 
SCDM model. Also, most clusters are closer to us in a SCDM model, thus the contribution 
from clusters to the power spectrum peaks at higher angular scales (lower £) than in low 
matter density models. The KSZ and MCG effects have their coherence length (peak 
contributions) at £ ^ 10000. The coherence length of the SSZ effect is about £ ps 1000 for 
our SCDM model and at about £ « 2000 for our OCDM and ACDM models. In general, 
the contributions to the power spectrum from the SSZ effect are about 2 and 3 orders of 
magnitude greater than those from the KSZ and MCG effects. At very small angular scales, 
£ > 10000, the contribution to the power spectrum from the MCG effect might exceed that 
of the SSZ or KSZ effects, and even the primordial fluctuations in the CMBR, but this 
depends on the details of the evolution of cluster atmospheres. 

Our simulations give somewhat different results for the KSZ and MCG effects than 
those of Aghanim et al. (1998). In our simulations the amplitudes of the KSZ and MCG 
effects for low matter density models are about an order of magnitude greater than those 
for our SCDM model, and have a coherence length of about £ = 10000, while rising 
monotonically at smaller I. According to Aghanim et al.'s simulations, with similar cut 
off to ours, Mi ow = 1O 13 M , contributions from the KSZ effect in all models constantly 
grow and show no signs of leveling off, and their amplitude has a very weak dependence on 
cosmological model. Contributions from the MCG effect on the other hand show a plateau 
in all Aghanim et al.'s models for £ > 1000, and for the SCDM model, the MCG effect is 
larger than for the other two models. 

Quantitatively, our models show cluster-related effects that are weaker by a factor of 
10 for the MCG effect in our SCDM model and a factor of 100 for the KSZ effect for all 
models. We attribute these differences mostly to the different evolution models for the 
intracluster gas. The ratio between the overall amplitudes of the KSZ and MCG effects in 
our calculations is about the same as in Aghanim et al. (1998) 's results. Our results show 
that the power spectra of the KSZ and MCG effects do not exceed the gravitationally lensed 
primordial power spectrum up to £ ~ 10000, while the power spectrum of the SSZ effect 
becomes dominant at £ ^ 7000 in all our models. Contributions to the power spectrum 
from the SSZ and KSZ effects based on Aghanim et al.'s model would exceed those from 
the CMBR at £ ^ 5000 even if one takes gravitational lensing of the primordial CMBR into 
account. Our results are comparable to those obtained by Tuluie et al. (1996) and Seljak 
(1996). Persi et al. (1995) 's results for contributions to the power spectrum from the SSZ 
(KSZ) effect are about the same as (an order of magnitude higher than) our results suggest. 
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8. Discussion 

An observed power spectrum is made up from the sum of all astrophysical effects and 
noise. We rely on the different frequency and/or power spectra of the secondary effects to 
separate these foregrounds from the primordial CMBR signal ( Tcgmark 1998 ). 



Of the effects discussed here, it should be easy to separate the SSZ effect by using 
multi-frequency measurements of its unique spectrum. The separation of primordial 
fluctuations in the CMBR and fluctuations caused by the KSZ and MCG effects is more 
difficult since their frequency spectra are the same. Optimal filters have been designed to 
separate the KSZ effect ( [Haehneft fc Tegmark 1996| ; [Aghanim et al. 1997| ): here it helps to 
know the SSZ effect for the same cluster, since that would give us a position and even an 
estimate for the expected amplitude of the effect. Aghanim et al. (1998) discussed methods 
to separate the MCG effect: this is facilitated by its unique dipole pattern with sharp 
peaks (Figure [|). Primordial fluctuations are usually assumed to be gaussian, where the 
probability of getting such a strongly peaked bipolar pattern is small, and we would expect 
the strong small angular scale gradient near a known cluster of galaxies to be a definite 
indication of the presence of the MCG effect. Also, knowing the position of clusters helps 
to find the effect. However, contributions from other effects, such as early ionization and 
discrete radio sources causes further confusion, and may be expected to make it difficult to 
determine the power from the SZ or MCG effects. 

We analyzed the contributions to the power spectrum from the SSZ, KSZ and MCG 
effects to check their impact on the determination of cosmological parameters, especially at 
large I where gravitational lensing may break the geometric degeneracy. In Figure |^ we show 
the small scale lensed primordial fluctuation power spectra of our three models (OCDM, 
SCDM, ACDM; solid lines) with power spectra resulting from the sum of fluctuations 
due to the lensed primordial CMBR and the SSZ effect (long dashed lines), and from the 
sum of the lensed primordial CMBR, the KSZ and the MCG effects (short dashed lines). 
According to our models, if the fluctuations due to the SSZ effect are fully separated, 
the KSZ and MCG effects do not prevent the use of this part of the power spectrum to 
break the geometric degeneracy and distinguish between different CDM models. Note that 
normalization at the first Doppler peak, rather than the usual COBE normalization, would 
lower the contributions to the power spectrum from primordial fluctuations in a ACDM 
model relative to those from a SCDM model, and thus secondary effects would become 
more important relative to the primordial CMBR fluctuations. 

Our simulations also show that the power spectrum of the SSZ effect may itself be 
used to break the geometric degeneracy. Since the separation of the SSZ effect from other 
secondary effects should be straightforward, we should be able to determine the power 
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spectrum of the SSZ effect alone. As can be seen from Figure |7], this power spectrum 
depends on Q h 2 and flbh 2 , providing an additional constraint on these parameters. 
We note, however, that the amplitude of the SSZ effect is model dependent. Since the 
contributions to the power spectrum from the SZ and MCG effects are model dependent, 
to evaluate fully their power spectra we need a better observationally-supported model 
for the intracluster gas. Sensitive, high-resolution all-sky, X-ray observations could map 
the emission from intracluster gas up to high redshift providing strong constraints on gas 
formation and evolution and thus a good basis for modeling the SSZ and KSZ effects 
flJahoda 1997] ) . Number counts of clusters based on their SSZ effect can also be used to 
constrain cosmological models QSuto et al. 1999| ). 

There are many possibilities of using observations to break the geometric degeneracy. 
For example measurements of the CMBR polarization, the Hubble constant, or light curves 
of Type la supernovae have been discussed (Zaldarriaga et al. 1997; Eiscnstcin, Hu fc 



Tcgmark 1998| ; |Tegmark et al. 1998Q . Also, combination of measurements of the SSZ effect 
and thermal bremsstrahlung (X-ray) emission from clusters can be used to determine the 
Hubble constant for a large number of clusters, providing a statistical sample which might 
enable us to determine the Hubble constant, and perhaps the acceleration parameter, to 
good accuracy ( [Birkinshaw 1998| ). 

Secondary fluctuations introduce non-gaussianity into the primordial spectrum at 
small scales. This non-gaussianity should be taken into account when estimating CMBR 
non-gaussianity at these scales. Winitzky (1998) estimated the effect of lensing and 
concluded that Planck may observe non-gaussianity due to lensing near the angular scale 
of maximum effect, ~ 10'. Other processes, including the SSZ, KSZ, and especially the 
MCG effect, introduce a highly non-gaussian signal as is easily seen for the MCG effect 
on Figure |3|. A similar non-gaussian pattern arises from moving cosmic strings (the 
Kaiser-Stebbins effect, compare our Figure |3| to Figure 6a of [Magueijo fc Lewin 1997] ) . 
Our results indicate that at £ ^ 10 4 the MCG effect might be comparable in strength to 
the primordial fluctuations. Evidence for non-gaussianity has been reported by Ferreira, 
Magueijo & Gorski (1998) and Gaztanaga, Fosalba, & Elizalde (1997) at angular scales 
£ sa 16 and £ ~ 150. They do not exclude the possibility that this non-gaussianity has been 
introduced by foregrounds, but our results show that clusters can not introduce detectable 
non-gaussianity on such scales (Figure |7|). 

We convolved our theoretical results (Figures f| - 0) with the expected point spread 
functions (PSFs) of instruments on the MAP and Planck missions to estimate the level of 
the secondary fluctuations caused by clusters of galaxies on the observable power spectrum. 
The observed Gt values become 
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Cf a = C e W e , (72) 

where the We values are the Legendre coefficients of the PSF. For an assumed gaussian 
response function and in the small angle approximation, the Wt coefficients are 



W e = e-^ e+1 ^ 2 , (73) 

where a = h/2\fhii, and h is the FWHM of the beam (for a detailed description of window 
functions and We coefficients, see White & Srednicki 1995). The observed rms fluctuations 
then become 



IP 4- 1 

(AT/T ) r 2 _ = £ -rr-CWt- (74) 

t 

In general, contributions from unresolved cluster static effects add to provide a 
cumulative contribution to the CMBR power spectrum. Contributions from the KSZ and 
MCG effects from unresolved sources tend to cancel. In the case of the MCG effect this 
is because each unresolved source contribution would be zero owing to the dipole spatial 
pattern of the effect. For small-scale KSZ effects there are several sources in the field of view 
of the telescope, and different sources have positive or negative contributions depending 
on the sign of their line of sight peculiar velocity, and therefore they tend to cancel each 
other. The larger the beam size, the more effective is the cancellation of the MCG and KSZ 
effects. Note that the spatial extension of the MCG effect is much larger than that of the 
KSZ effect, so many clusters may be unresolved in their KSZ and resolved in their MCG 
effect. The MCG effect might be relatively more important at high redshifts, since it does 
not require a well-developed cluster atmosphere. 

In Figures [| and |10| we show the contributions to the power spectrum from primordial 
fluctuations, and the SSZ, KSZ and MCG effects, convolved with the PSF of the planned 
v = 94 GHz receiver on MAP, and the planned v = 353 GHz bolometer on Planck. The 
amplitude of the fluctuations from the the SSZ effect is negative at v = 94 GHz and 
positive at v — 353 GHz, but only the the absolute value of the effect contributes to the 
power spectrum. The different maximum t values for the MAP and Planck systems (£ ma x ~ 
1000 and 2000, respectively) can clearly be seen on Figures |9] and |H]. Because of these 
cutoffs, the observable power spectrum is dominated by primordial fluctuations at all £ for 
these missions. According to our results, the SSZ effect may cause a 1% enhancement in 
the amplitude of the Doppler peaks, which is at the limit of the sensitivity of the MAP 
and Planck missions. From the analysis of the power spectrum, this would lead to an 
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overestimation of the parameter Q h 2 by about 1%,. The shift in the position of peaks as a 
function of £ caused by the SSZ effect is less important since the spectrum of the SSZ effect 
has only a weak dependence on £. 

In Table [I] we show the (AT/T) rms values of the contributions to the CMBR from the 
SSZ, KSZ, and MCG effects convolved with the the 94 GHz MAP and 353 GHz Planck 
receivers for our three models (OCDM; ACDM; SCDM). As a comparison, we display 
the corresponding rms values of the primordial fluctuations. The rms values of all these 
secondary effects are an order of magnitude smaller than rms values from primordial 
fluctuations. The most important contribution is that of the SSZ effect at these frequencies. 
The KSZ and MCG effects give similar contributions with the KSZ effect being about a 
factor of two stronger. Aghanim et al. (1998)'s results for the rms values of the MCG effect 
is a factor of 10 (SCDM) or a factor of 3 (OCDM and ACDM) larger than our results. 
Note, however, that rms values give only a crude estimate of the magnitude of the effects. 
At large angular scales the primordial fluctuations are about 100 (for the SSZ effect) or 
10 5 — 10 7 (KSZ, MCG effects) times stronger than the secondary fluctuations. 

An ideal observation to measure the contribution to the power spectrum from the SSZ 
effect would use high angular resolution (£ ^ 7000) and high frequency [y ^ 250 GHz). 
SuZIE probes the power spectrum at angular scale i « 7500 at 140 GHz. The 2a upper 
limit on the power at this scale from SuZIE is £{£ + l)Cg < 1.4 x 10~ 9 flGanga et al. 1997f) . 
Unfortunately our models suggest that at this frequency the primordial contribution to the 
power spectrum is about 10 times stronger than that from the SSZ effect. A promising 
experiment is SCUBA, which probes the anisotropies at angular scale £ ~ 10000 and 
frequency 348.4 GHz. Their preliminary 2a upper limit on the power spectrum at this scale 
is £{£ + l)Cg < 4.7 x 10~ 8 . Much further work is planned, and should lower this limit by a 
factor of 3-10 ( [Borys, Chapman fc Scott 1998| ). 
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Appendix: Small Angle Approximation for P\ 

We derive a small angle approximation for the associated Legendre polynomial P}. We 
express the associated Legendre polynomial Pj in terms of the Legendre polynomial P^ (see 
for example |Arfken fc Weber 1995| ) as 

PKfi = V /rr ^^^(^)- (!) 

Introducing a new variable, x = £(£ + using the chain rule and a small angle 
approximation, we get 



Using a small angle approximation for Pi (Peebles 1980) we obtain 



-£(£ + l)^*-£(£+l)-±J ( x ). (3) 
ax ax 



Using the relation 



Ji(x) = -j^Mx) ( 4 ) 



for Bessel functions, we finally obtain 



P/(/x) « -£(£ + l)^Jo(x) = £(£ + 1) Ji(x). 



(5) 
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Fig. 1. — The PSMFs, dN / (dV d\n M)(> M), at z — 0, used in our simulations compared 
to observations. The long dashed, dash dot, and short dashed lines refer to the PSMFs of 
our A = 0, tt = 1 (SCDM); A = 0, fi = 0.2 (OCDM); and A = 0.8, fi = 0.2 (ACDM) 
models respectively. The solid line represents the observationally-derived mass function of 
Bahcall and Cen (1993). 
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Fig. 2.— High mass (M > 1.6 x 1O 15 M ) PSMFs, dN/(dVdz) in units of number Mpc~ 3 , 
as a function of redshift, z, used in our simulations compared to observations. The long 
dashed, dash dot, and short dashed lines refer to the cumulative PSMFs of our A = 0, 
tt = 1 (SCDM); A = 0, tt = 0.2 (OCDM); and A = 0.8, fi = 0.2 (ACDM) models 
respectively. The solid line represents the mass function derived from observations (Bahcall 
et al. 1997), which is a close match to our OCDM model. 
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Fig. 3. — Fluctuations in the CMBR (AT in fiK) caused by the moving clusters effect in our 
SCDM model (see text for details). Primordial fluctuations and other source of secondary 
fluctuations are not displayed. The area covers 30° x 25° of the sky with pixel size 2.6' x 2.6'. 
The maxima and minima are about ±2/xK. 




Fig. 4. — The power spectrum C(£) = 1(1 + l)Ce of the lensed primordial CMBR (solid 
line), SSZ (long dashed), KSZ (dashed dot), and MCG (short dashed) effects for an OCDM 
model with Q = 0.2, A = 0, n = -1.4. 
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Fig. 5. — The power spectrum C(£) = 1(1 + l)Ce of the lensed primordial CMBR (solid 
line), SSZ (long dashed), KSZ (dashed dot), and MCG (short dashed) effects for a ACDM 
model with Q = 0.2, A = 0.8, and n = -1.4. 




Fig. 6. — The power spectrum C(£) = 1(1 + l)Ce of the lensed primordial CMBR (solid 
line), SSZ (long dashed), KSZ (dashed dot), and MCG (short dashed) effects for an OCDM 
model withQ = 1, A = 0, and n = —1.4. 
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Fig. 7. — Power spectra C(£) = 1(1 + l)Cg: lensed primordial CMBR (solid line), SSZ (long 
dashed), KSZ (dashed dot), and MCG (short dashed) effects for our three models: OCDM, 
ACDM, and SCDM. 
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Fig. 8. — Small scale power spectrum C{£) = £{£ + 1)Q of the lensed primordial CMBR, 
the static and kinematic SZ effects, and the MCG effect for our three models: OCDM, 
ACDM, SCDM. The long and short dashed lines represent the sum of the lensed primordial 
CMBR and the SSZ effect, and the sum of the lensed primordial CMBR, the KSZ and RSC 
effects respectively. The contribution from the static SZ effect is important in all models. 
Contributions from the kinematic SZ and MCG effects are negligible at these I values. We 
may conclude that if fluctuations due to the static SZ effect are removed the kinematic and 
MCG effects do not prevent using gravitational lensing to break the geometric degeneracy. 
Note that different normalizations shift the sets of curves. For example, normalization to the 
first Doppler peak makes the power spectra of our SCDM and ACDM models shift closer. 
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Fig. 9.— Power spectra C(£) = £(£ + l)C t of fluctuations in the CMBR as MAP would 
observe them at v = 94 GHz with FWHM = 12.6'. Models and line codes are the same as 
in Figure 0. The contribution from the primordial fluctuations dominates on all scales. 
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Fig. 10. — Power spectra C(£) = £(£ + l)Ce of fluctuations in the CMBR as Planck would 
observe them at v — 353 GHz with FWHM = 5 arc minute. Models and line codes are 
the same as in Figure [7[ The contribution to the power spectrum from the static SZ effect 
at this frequency is 1.2 times that in the Rayleigh- Jeans region, and the static SZ effect 
dominates over the processed primordial fluctuations above about I ~ 3000 (except for the 
SCDM model). 
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Model 




MAP 


94 GHz 






Planck 


353 GHz 






CMBR 


ssz 


PKSZ 


MCG 


CMBR 


SSZ 


PKSZ 


MCG 




(jtK 




(JtK 


liK 


fiK 


/zK 


/j,K 


fiK 


1. OCDM 


86 


7.0 


0.54 


0.20 


100 


14 


0.81 


0.33 


2. LCDM 


119 


5.9 


0.49 


0.20 


130 


11 


0.71 


0.32 


3. SCDM 


93 


3.2 


0.36 


0.12 


101 


5.7 


0.50 


0.17 



Table 1: Root mean square (AT rms ) of the SSZ, KSZ, and MCG effects convolved with 
response functions of the 94 GHz MAP (FWHM = 12.6') and 353 GHz Planck (FWHM 
= 5') instruments from our three models (OCDM; LCDM; SCDM, see section [l]). As a 
comparison, we display the rms values of the primordial fluctuations (CMBR). 



